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We derive an extended lattice gauge theory type action for quantum dimer models and relate it 
to the height representations of these systems. We examine the system in two and three dimensions 
and analyze the phase structure in terms of effective theories and duality arguments. For the 
two-dimensional case we derive the effective potential both at zero and finite temperature. The 
zero-temperature theory at the Rokhsar-Kivelson (RK) point has a critical point related to the 
self-dual point of a class of models in the N 00 limit. Two phase transitions featuring a 
fixed line are shown to appear in the phase diagram, one at zero temperature and at the RK point 
and another one at finite temperature above the RK point. The latter will be shown to correspond 
to a Kosterlitz-Thouless (KT) phase transition, while the former will be governed by a KT-like 
universality class, i.e., sharing many features with a KT transition but actually corresponding to a 
different universality class. On the other hand, we show that at the RK point no phase transition 
happens at finite temperature. For the three-dimensional case we derive the corresponding dual 
gauge theory model at the RK point. We show in this case that at zero temperature a first-order 
phase transition occurs, while at finite temperatures both first- and second-order phase transitions 
are possible, depending on the relative values of the couplings involved. 
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I. INTRODUCTION 

The "quantum dimer model" (QDM)ii^ was intro- 
duced to emulate the quintessential features of valence 
bond states in spin systems. Its inception was moti- 
vated by the short-range resonating valence-bond (RVB) 
state of Andersouj^ which was used as a possible start- 
ing point for a theory of high- Tc superconductors 4 Since, 
the QDM has played a prominent role in modeling vari- 
ous frustrated magnets,^ cold atom systems,^ Josephson 
junction arrays;^ spin-orbital systems;^ topological quan- 
tum orders?^ fractionalization, and other phenomenal^ 
The only degrees of freedom in the QDM are dimers that 
represent singlet states formed by neighboring spins on 
the lattice. Two spin singlet states cannot overlap or 
share a common lattice site. Similarly, any two dimers 
within the QDM cannot overlap — the dimers satisfy a 
"hard core" constraint. Several spin^ i^°'^^i^^ , and orbital 
models^ indeed have ground states which are precisely of 
the dimer type. In several spin systems, pairs of 5* = 1/2 
spins bind into singlet states. In some spin systems, such 
as the Klein models,-*^ it can be proven that not only 
are dimer states ground states but that they are the only 
ground states on general lattices .J^ 

One of the major complications of real spin sin- 
glet systems by comparison to the QDM, is the non- 
orthonormality of the singlet product basis states (an 
item which should not be confused with viable linear in- 
dependence of these states). In spin systems with dimer 
ground states is possible to systematically write down 
rules for the evaluation of the overlap between the sin- 
glet states on bipartiteii^ii^ and more general lattices^i 
as well as to evaluate the general matrix elements of spin 



exchange and other terms.— The QDM avoids many 
complications by focusing on the quintessential physics of 
hard core dimer systems. The overlap between different 
dimer states is simply set to zero. 

The Hamiltonian of the QDM on a square lattice 
readsi 

H = EH(III)("I + I«)(III) 

□ 

with the sum performed over all elementary plaquettes of 
the square lattice. We reiterate that the dimer states in 
this model are orthonormal,^ differently from the singlet 
valence-bond states originary from quantum spin models, 
which are in general not orthogonali^ 

As seen in Eq.(IT]), the QDM contains both a kinetic 
(t) term that flips one dimer tiling of any plaquette to 
another (a horizontal covering to a vertical one and vice 
versa) and a potential (v) term counting how many pla- 
quettes are flippable. These two terms — kinetic and 
potential — were treated on different footing in earlier 
works. Particularly important is the point t — v, the 
so-called Rokhsar-Kivelson (RK) point, where both the 
kinetic and potential terms are of equal magnitude. The 
RK point is an exactly solvable critical point that sepa- 
rates two different valence-bond solid (VBS) phases. 

Within the VBS phases, the dimers break lattice point 
group symmetries. ^'^^ On non-bipartite lattices and t ^ 
V, valence bond liquid phases appear. On the triangu- 
lar lattice, short ranged resonating valence bond phase 
with no gapless excitations and with dcconfined, gapped, 
spinous appear for a finite range of parameters. Similar 
behavior is also found on the kagome lattice. For square 
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lattices, valence-bond liquid phases seem to be possible 
only at the RK point. 

It has been shown in several recent 
paper o^'''^°i^^'^^i^^'^'^ that the phase structure of 
the QDM is considerably richer than previously thought. 
Especially interesting here is the emergence and char- 
acterization of the different VBS phases. In particular, 
it has been recently shown^^ that the VBS phase 
continuously interpolates with the plaquette phase via 
a mixed regimei^ Classical three-dimensional dimer 
models have also recently shown to exhibit interest- 
ing phase transitions. Similar to some quantum 
phase transitions in antiferromagnets^^^ several three- 
dimensional classical dimer models seem to exhibit 
second-order phase transitions that do not easily fit 
within a Landau- Ginzburg- Wilson approach.^'* 

In this paper we will study the QDM of Eq. ([T]) via ef- 
fective field theories and duality methods. Our work will 
focus on the square and cubic lattices. An important goal 
of this work is to elucidate the nature of the phase transi- 
tion, both at zero and at finite temperatures. We charac- 
terize the phase transitions in two and three dimensions 
by employing an interplay between renormalization and 
duality. The plan of the paper is as follows. In SectionlTTl 
we will introduce a simple spin 5=1/2 representation of 
the QDM which will allow us in Section lTTTl to introduce a 
new lattice gauge theory representation of the QDM. We 
show that the QDM is equivalent to an extended Abelian 
lattice gauge theory in which in addition to the usual 
plaquette cosine term of standard lattice gauge theories, 
a higher harmonic of the field strength is also present. 
It is known from the lattice gauge theory literature^^ 
that in four spacetime dimensions such a theory has a 
rich phase structure. In Section HVl the two-dimensional 
QDM model is discussed using an effective height model. 
Such height model representations of the QDM are well 
known in the literature,— '^Si'^'^i^ and have been moti- 
vated by physical arguments. Here, we put these repre- 
sentations in perspective via the extended lattice gauge 
theory that we derive in Section IIIII We show that the 
height field theory model can be derived from the ex- 
tended lattice gauge theory. It should be stressed here 
that the extended lattice gauge theory in Section II is 
more precise than the height models discussed in Section 
IIVI However, studying it by purely analytical means is 
very difficult, so that effective height models are indeed 
useful in this respect. In Section |Vl the phase structure 
of the two dimensional QDM is analyzed using effective 
potentials and duality arguments. The phase diagram 
is discussed both at zero and finite temperatures. We 
show that the full phase diagram features two Kosterlitz- 
Thouless (KT) like phase transitions, one at zero temper- 
ature at the RK point, and another one at finite temper- 
ature away from the RK point at t > v. The latter is a 
genuine KT transition, while the former is KT-like, i.e., 
it shares many properties with the usual KT transition, 
but it actually corresponds to a new universality class 
featuring a fixed line. We further illustrate that the ef- 



fective potentials may be the same in both situations. In 
Section rvn we discuss the three-dimensional QDM at the 
RK point at both zero and finite temperatures. To this 
end, the extended lattice gauge theory is considered more 
directly in a Villain approximation. Here duality plays 
a crucial role in determining the phase structure of the 
RK point. We show that at zero temperature and at the 
RK point a first-order phase transition between a VBS 
and a RVB state takes place. Thus, in contrast with the 
two-dimensional QDM, no quantum critical point exists 
in three dimensions at the RK point. For finite temper- 
atures both first- and second-order phase transitions are 
possible. The character of the transition (whether it is 
continuous or abrupt) depends on the relative values of 
the couplings involved. We summarize our conclusions in 
Section Eni 

II. SPIN REPRESENTATION OF THE QDM 

In this section, we will find a direct representation for 
the QDM in terms of a lattice gauge theory. This will al- 
low us to systematically derive the height representation 
as an approximation. To achieve this aim, we first express 
the QDM directly in a spin language. As we will show 
in a future publication, the approach that we introduce 
below enables a derivation of spin representations and en- 
suing gauge type theories for other lattices with general 
(non-square type) elementary plaquettes. We will employ 
a simple algebraic property of the QDM: The potential 
term in Eq. {Ip is the square of the kinetic term, 

+ (2) 

What enables this relation is the orthonormality of the 
dimer states.»i The relation of Eq.(l2]) will enable us to 
treat both the kinetic and potential terms within a uni- 
form systematic framework later on. All (positive) even 
powers of the kinetic term give rise to the potential en- 
ergy term whereas all odd powers of the kinetic energy 
term yield the kinetic term unchanged. Denoting the 
kinetic term by B\j, Eq.([T]) can be rewritten as 

i7 = ^(-ii?n + i;S^), (3) 
□ 

Thus, for a given value of w > 0, the ground states min- 
imize the sum — t/i'^v))'^ over all plaquettes of 
the latttice. As will become evident later on, the kinetic 
term _Bq will play, in the gauge representation, a role sim- 
ilar to that of a (modular) magnetic flux that threads a 
plaquette. The flux pattern Bu within the ground states 
will be determined by the ratio t/v. 

We will now employ and cast the relation of E c|.(|p in 
a spin language. Similar to Refs. O, HI, and |34 we 
designate the presence/absence of a dimer between the 
two sites i and j by erf = 1 and erf - = — 1 respectively. 
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Introducing the Pauh raising/lowering operators 
i {(jfj ± icr^j ) , the QDM Hamiltonian reads 



H = ~tJ2iWa + W^) + vY,iWaW^ + W^Wn), (4) 



where 



(5) 



and □ = ijkl is a plaquette. The spin representation 
of the kinetic only term was discussed in Refs. [29ll33ll3^ . 
Our new full Hamiltonian of Eq. containing both the 
kinetic and potential terms, is a faithful representation 
of the QDM of Eq. (P). It automatically incorporates 
the hard-core constraint of the QDM that prevents two 
dimers from overlapping. If we denote the two verti- 
cal bonds of the plaquette |J J), by (ij) and (kl), then 

1 1 I) — Im)- ^'^^ state orthogonal to |j J), 
we have W^\ip) = 0. Similarly, W^n IHI) = ll I)- The 
products (W^nW^g) and (VFgVFn) constitute projection 
operators onto the states 1 1 J ) and | ^ ) respectively. As 
(CTj^)^ = (cTjj)^ = 0, we obtain the Hamiltonian of Eq.Q 
where 



(6) 



The kinetic term in Eq. (j4|) is precisely the ring exchange 
term in a well studied model (in its XY version)35- in 
which we regard the centers of bonds of a square lattice 
as vertices of a square lattice rotated by 45 degrees and a 
scaled down lattice constant by a factor of 2~^/^. In the 
large spin limit, Eq.Q is a classical XY Hamiltonian. By 
introducing a gauge field representation of the spins, we 
can check whether deconfined criticality^^ may arise. 



III. LATTICE GAUGE THEORY 
REPRESENTATION 



We will now transform the exact spin representation 
into a lattice gauge theory by writing down the spin co- 
herent basis action corresponding to the Hamiltonian of 
Eq. ([3]). To this end, the spin (s — ^a) points anywhere 
on a sphere of radius s — 1/2. In imaginary time, the 
spin performs a cyclic evolution on the spin during the 
time interval < r < /3. The Euclidean action is 



Parameterizing the spin on the bond {ij) 
by its location on the sphere by Sy = 
s{sm9ij cos Aij , sin 9ij sin Aij , cos 9ij), the Hamilto- 
nian of Eq. ^ reads 

H = 2J2 [-tisind)a cos(F™„) 

■i,m,n 

-ht;(sin6l)^cos(2F,™„)], (8) 

where Fimn = '^m.Ain — Vn-A-im is the field strength in 
the lattice, with the latin indices m, n being summed over 
(d — 1) spatial dimensions, and i a site label. The fac- 
tor (sin0)n denotes the product of four factors of sin6'ij 
around a plaquette. The area subtended by the spin on 
the sphere depends on Oij. The WZW contribution can 
be rewritten as 



I cost 



dr 



(9) 



The height representation of the QDM was introduced 
and motivated by numerous insightful intuitive consider- 
ations in earlier work o^'^^i'^^'^^i'^'^'^^ . In the present work, 
we derive it by a sequence of approximations. 
The partition function 



Z 



(10) 



with S the action of Eq.Q. Henceforth, we will make 
several approximations. In the Hamiltonian, we will re- 
place the explicit angular dependence on the angles 9ij 
(to be integrated over) with an average value determined 
(within the coherent path integral formulation) by set- 
ting \sfj\ = 1/2. Later on, we will further invoke both 
the continuum and Villain approximations. 

Within the coherent path integral formulation, s^^ is 
continuous. We approximate the integral over the po- 
lar angles 9ij by enforcing the S = 1/2 results that 
would have been obtained in the canonical formulation 
(i.e. that would have been obtained for the averages of 
these squares alone with the coherent spin path integral) . 
To this end, we replace in the Hamiltonian any appear- 
ance of sin^ 9ij by 1/2, as the squared norm of the XY 
(planar) part of the S = 1/2 spin is s{s + l)—{s^)^ — 1/2. 
That is, in all pertinent expressions for the spin 5* = 1/2 
system we set |sf | = 1/2 and 



S = -is ^(2^ - A^j) + I dr H[{c7ij}], 

ij 



(7) 



where in the first (Wess-Zumino-Witten (WZW) borne) 
term, Aij is the area subtended by the spin Sij — ^aij 
as it cyclically evolves in imaginary time (0 < r < /3) on 
the sphere. This term corresponds to the Berry phase 
within the coherent spin basis. In Eq.([7]), {s^} denotes 
all spins of the lattice (all bond centers). 



(11) 



with real compact . Substituting this form in Eq. (O, 
we obtain, from Eq. ([3]) , the Hamiltonian of an extended 
lattice gauge theory, 



t 



cos{Fi, 



,) + - cos{2Fi, 



(12) 
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We can construct an action in which crf^ plays the role 



Lagrangian density of the form, 



of a "momentum" that is conjugate to the phase Ai 



S — is a 



\7 A - 



H. 



(13) 



In Ea. (|13p . we discretize the imaginary time r. The first 
term with s(t|„ — > s|„ is the WZW term of Eq. ^ ~ 
it is determined by the area A swept by the spin Sj„ as 
it traverses the sphere in the closed orbit for < r < 
[see Eq. ©]. This Berry phase is similar to the one 
introduced in the Sachdev-Jalabert model of quantum 
antiferromagnetsi^ The WZW term may also be red- 
erived from another vista. As [crf^ , c,^] = ='=2o'j"^, it fol- 
lows that the lattice gauge field Aij is canonically conju- 
gate to af^ (this is similar to the discussion in Ref. 38: in 
a different context). The partition function is therefore 
given by 



Z = J Yl '^^i" *5^P m '^jn'^rAjn - H 

in {a' =±1} \ T,j,7i 

(14) 

As stated earlier, within the coherent path integral for- 



mulation, 



is continuous. To be consistent, we in- 



voke the same approximation that we employed earlier 
{\sfj\ = 1/2) and set cr^ = ±1. Summing over afj, this 
leads, up to an irrelevant constant, to the result 

Z = j WdA,n Y[ cos{sVrAjn)exp{-H) 

in T-j-jTl 

[]dA,„e-^, (15) 



where the action in terms of the gauge field only is given 
by 



^ = - ^ ln[cos(sV^A,„)] 



i.T.rn.n 



L o 



(16) 



A related theory in which the first term in Eq. ([T| 
was absent was studied numerically in four dimensions.^® 
Such a theory usually has a rich phase structure char- 
acterized by a tricritical point. A classical XY version 
was also studied in three dimensionsi^ Note that the ex- 
tended lattice gauge action of Eq. (|16p differs from earlier 
proposals for the QDM.^'''''* Amongst other benefits, our 
lattice gauge action enables us to derive the height model 
as an approximation without resorting to phcnomcnolog- 
ical arguments. It also allows us to immediately identify 
t = u as a special point. 

Expanding the action (jl6p up to quadratic order in the 
fields, taking the continuum limit, and rescaling both i 
and by a uniform factor, t — s- s^t, v — > s^i;, we ob- 
tain up to an overall innocuous multiplicative factor, a 



t - V 



(17) 



When we use our quadratic gauge action, we have a de- 
pendnce only on (i — v) and not on the individual values 
of t and V. This symmetry in the quadratic order gauge 
(and, as we will see later on, the derived height model) 
description of the QDM is lifted as higher order terms 
are included in our gauge action. 

By analytically continuing to real time, we have the 
following gauge field propagator in the transverse gauge 
V • A = 0, 



£'„„(cj,p) = 



1 



uj"^ — [t — i;)p2 



c PmPn 



(18) 



We see that ioi t > v the spectrum of elementary excita- 
tions is given by 



Vt^\p\. 



(19) 



For t < V, on the other hand, we have a purely damped 
mode. This is not a completely satisfactory physical re- 
sult, since it makes the system unstable, especially at 
finite temperature. In order to stabilize the system in 
this case, we have to include higher order terms, which 
can be obtained from the photon self-energy for uj small. 



IV. DERIVATION OF THE 
TWO-DIMENSIONAL HEIGHT MODEL 

The quintessential low-energy properties of the QDM 
can be captured by a height representation2i2£i2ii2^iMi2& 
where the weights (heights) h are associated with the 
plaquettes of the underlying lattice. Our gauge theory 
action of Eq. (|16p — an average coarse grained form of 
the exact QDM Hamiltonian in the form of Eq. ^ — 
offers a new route towards systematically obtaining the 
height model representations of the QDM. In this sec- 
tion, we will show to obtain these representations from 
the gauge theory of the previous section. The advan- 
tage of the lattice gauge theory approach is that it al- 
lows to easily obtain the field theories in both the two- 
and three-dimensional cases. In a nutshell, the gauge 
theory and height representation are related by a simple 
( "Abelian-Higgs" type) duality. This duality amounts to 
the following rule of thumb in 2+1 dimensions. 



ii^xd^Ax oc df_,h, 



(20) 



with h the height field and where the greek indices run 
over the spacetime coordinates. That is, an effective mag- 
netic field derived from our gauge potential is equal to the 
gradients of the height field. Note that due to the com- 
pactness of the gauge field, the magnetic fiux through a 
closed surface should be quantized. 

To make it more lucid, we will briefly describe the 
physical content of the correspondence of Eq. (|20p for 
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the cases represented in Fig. [2l We simply focus on the 
relation of Eq. ([^0]) to see to what height phases the 
various VBS phases correspond to. First, we define the 
flux per plaquette $ as the counter-clockwise sum of the 
(directed) link variables Aij . 

(a) Inspecting Eq. we see that when $ = (mod 
2tt), no flux pierces the plaquettes and by Eq. (|^D|) . the 
height field /i is a constant. This case corresponds to 
the so called plaquette phase, where the dimers resonate 
around a plaquette in a would-be columnar pattern. 

(b) Similarly, the staggered VBS corresponds to half a 
flux quantum per plaquette 



(mod 2tt). 



(21) 



With a fixed circulation direction fiux in Eq. (jSTj), this 
flux is seen to be staggered on the two sublattices for a 
flxed value of Aij . In the case of Aij = 7r(mod 2it) as tt = 
— 7r(mod 2tt), a staggering of half-fluxons is equivalent to 
a uniform array of half-fluxons (sans staggering). 

Both states (a) and (b) are invariant under the time 
reversal operation which corresponds to A^j — > —Aij. 

To derive the correspondence of Eq. we will follow 
a sequence of steps below, (i) We will write down a Vil- 
lain type action that reduces to Eq. (|17p in a continuum 
limit, (ii) We will then apply the Poisson summation 
twice on this action, (iii) Finally, we will show that the 
continuum limit of (ii) leads to a height model with the 
substitution of Eq. f^ . 

(i) We start by introducing the following effective lat- 
tice action in the Villain approximation: 



Like all Villain actions, it has the advantage of being 
more tractable, and yet including all relevant physics of 
the problem. 

Note that by considering a naive continuum limit and 
low momentum such that D(xi — Xj) « Sij/ p and rescal- 
ing the gauge field as Am \fP-^rm with p cx (t — u), 
we obtain a Lagrangian of the form (|17p . This low- 
momentum behavior reflects the fact that neglecting the 
compactness of the gauge field in Eq. (122]) , a Bogoliubov- 
like spectrum is obtained. 



c cK 



(25) 



For p > and p small we approximately recover the 
spectrum ([TO)) . We see that for p < a momentum space 
instability occurs. We will come back to this point later. 

(ii) We next dualize the action of Eq. ([22|) . By using 
the Poisson summation formula we rewrite the action in 
terms of auxiliary integer fields MjQ and Mj^ in the form 



i 

— i 'y^X-^'IjO^Omn'^mAjn + MjratmOn VrA,r^)■{2Q) 



We now use partial summation and integrate out the 
compact gauge field to obtain the constraint 



S = i^I?(x,-x,) 

X (Vt-Ajto — 2t: LijTn){SJ T A jjj^ — 27rZ/j,„) 



where 



The constraint is solved by 

Mjo = VrNj, Mj = VN. 



(27) 



(28) 



where Nj is an integer field. This leads to the action 



Z?(x, - X,) = 1 ^ e'P(''--^)7?(p), (23) S=\Y. 



(29) 



with 



^(P) 



[4 — 2(cospi -I- C0SP2)] + p' 



(24) 



where Lim, Ni^n 



are integer fields and p oc t ~ v. The 
parameter c which will appear in the height model that 
we will derive is usually set to unity in the literature. 
However, this should be avoided, as c is not dimension- 
less. In fact, as we will see, it plays a crucial role in the 
characterization of the quantum phase transition of the 
theory. 

The Villain gauge theory action constitutes a ver- 
sion of the original lattice gauge theory (|16p of the QDM. 



By using the Poisson formula once more to convert the 
integer field Nj into a real field hj , we obtain 



S = 



c{Wrh,f+p{Wh,f + ^{W'h,) 



27ri rijhj] 



(30) 



where Uj is a new integer field. 

(iii) In the continuum limit of the action (j30p the 
"charges" rij = ±1 are the most relevant ones, so that in 
the grand-canonical partition function the (product in j) 
of factors exp(i2TTnjhj) exactly exponentiates,^*' produc- 
ing a term cos(27r/i) in the continuum limit of the action. 
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TABLE I: Summary of the phase structure of the two- 
dimensional height representation of the QDM for different 
temperatures T and p oc (t — v). Kosterlitz-Thouless 
(KT) type transitions occur when (i) T = p = or (ii) 
T > 0, p > 0. At high temperatures, the system exhibits 
algebraic correlations for all p. 





p = 


p > 


p < 


T = 


KT-like 


Plaquette, VBS 


Staggered VBS 


T > 


No transition 


KT 


VBS melting 



Thus, we now finally arrive at a directly derived form the 
imaginary time Lagrangian of the height model 



C = ^{drhf + ^{\/hf + —{y^hf ~ZC0s{2TTh). (31) 
Z Z ZK 



Comparing Eg. ipT]) with Eq. (P^ . the correspondence 
of Eg. ((20)) becomes manifest. The height model of 
Eg. pip has been introduced and motivated by nu- 
merous insightful intuitive considerations in earlier 
works i^i^i^i^iMil^ Our approach enables a direct deriva- 
tion of the height model from the original QDM Hamil- 
tonian of Eq.([T]). 

Putting all of the pieces together, in 2-1-1 dimen- 
sions, our lattice gauge theory can be rewritten as a 
Lagrangian of Eg. pTj) for a scalar field, with the term 
— zcos(27r/i) following from the compactness of the lat- 
tice gauge field. 

In the case of 3+1 dimensions which we will dwell on 
later, the dual field strength is a second rank tensor, just 
as the original field strength. Therefore, we cannot intro- 
duce a scalar field in this case and have to work further 
with a gauge field. We will consider this case later on in 
detail (Sect. IV). 

The action of Eg. ((30| may, in principle, also lead to 
higher cosine harmonics, like for example cos(27rm/i) with 
TO — 2, 3, ... . However, as long as the couplings in 
front of all factors — cos(27rTO/i) are positive, the higher 
harmonics are irrelevant in the Renormalization Group 
sense.— On the other hand, there are situations where 
couplings with a negative sign play an important role, 
like in the case of emerging plaquette and mixed phases, 
as discussed recently in Ref. \2% 



V. TRANSITIONS IN THE 
TWO-DIMENSIONAL SYSTEM 

In this section, we study the structure of the phase 
diagram of the height model of Eq. ((3T|) . The outcome 
of our analysis is summarized in Table I. 

In what follows, we will examine the height model rep- 
resentation of Eq. (pij) for different values of p cx (t — u) 
and temperatures T. In terms of the effective theory of 
Eg. ipTj) . the RK point corresponds to p = 0. For p > 



the gradient term (V/i)^ dominates over (V^/i)^ at large 
distances, so that the latter can be neglected. The free 
propagator of the height field in momentum space and 
real time is 



Go{uj,p) = {h{uj,p)h{~uj,-p))Q 
1 



PP 



K 



(32) 



which diverges for an excitation spectrum of the form 
(l25)l . As the large distance limit is equivalent to p 0, 
we see that for p > the p'^ term is negligible. By in- 
tegrating out hj in Eq. ([50]) . we obtain a charged gas 
with charges Qj — 2Tmj interacting via a potential whose 
continuum limit yields the Euclidean counterpart of Eq. 
([32]) . For if = oo, in which case the term proportional to 
is also absent, this is just a classical Coulomb gas in 
three dimensions. It is well known that in this case the 
excitations are gapped )^^'^^ with a gap that never van- 
ishes. This result remains valid for K finite. Invoking a 
simple mean- field screening argument, the Debye-Hiickel 
approximation, and truncating the cosine term in (j3ip to 
lowest (quadratic) order, we arrive at the propagator 

G{uj,p) = {h{uJ,p)h{-LJ,~p)) = —3- — -. 

Gq {uj,p)-Att^z 

(33) 

It is readily seen that this propagator screens the inter- 
action between the charges of the gas. We turn back and 
ask what phases are to be found at different values of p. 

Essentially there are two types of phases. There are 
valence-bond phases where the height field is uniform and 
others where it modulates with some wave vector p^. The 
former leads to a height susceptibility that diverges for 
Lo = Q and p = 0, while the latter diverges for p^ = 
Pq when the frequency vanishes. In detail, we have the 
following scenarios: 

(a) For p > we can set if = cx3, since in this case 
higher gradient terms are unimportant. In this case we 
have to distinguish between the cases z < and z > 0. 
For z > 0, the height susceptibility ([33]) diverges only 
for Lu = and p = 0. This corresponds to the plaquette 
state.42, For z < 0, on the other hand, there is a diver- 
gence for a nonzero wavenumber po given by 

— • (34) 



This state corresponds to a VBS. These two phases — the 
plaquette and the VBS (flat) phases — are schematically 
shown as a function of p in Fig[T] Note that the plaquette 
phase occurs first, after the RK p oint (p = 0), before the 
VBS phase; see Refs. [13 andl25l 

(b) For p < the higher gradient term becomes impor- 
tant at large distances as a way to stabilize the system. 
In this case, the zero-frequency height susceptibility di- 
verges for a nonzero wavenumber po given by 



(35) 
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This leads to a modulation in the height field associated 
to a staggered VBS state. Note that for K — oo Eq. ([55]) 
reduces to Eq. ([34|) . with the difference that now it is p, 
and not z, that is negative. 

Below we will describe different regimes in the QDM 
both at zero and finite temperatures. In what follows it 
will be more convenient to discuss the finite temperature 
case first, since this essentially corresponds to a classical 
dimer model. 



A. Finite temperature, p > 

For p > and high temperatures (T), we effectively 
have a classical dimer problem, as in this case /i(r, x) « 
/i(x) and the system becomes effectively two-dimensional. 



S 



dr I d?xC « T-^ / d?xC 



d xL. 



classical . 



(36) 

After some trivial rescaling, we obtain the classical dimer 
Lagrangian 



-^classical ^ (^^) 



C cos(27r^/K/i), 



(37) 



with K — Tj p and ^ — z/T. In this limit, we ob- 
tain an ordinary two dimensional sine-Gordon theory4i 
This theory describes the vortex unbinding transition in 
the Kosterlitz-Thouless (KT) phase transition.'*^ It fol- 
lows that, the high-temperature regime of the RK model 
for yO > undergoes a KT phase transition^ in which a 
VBS or a plaquette stated melts into a liquid of valence- 
bonds, i.e., a finite temperature "spin liquid". The La- 
grangian (j37p corresponds to the finite temperature dual 
model^ to compact Maxwell electrodynamics in 2 + 1 
dimensionsi^ Although at zero temperature no phase 
transition occurs in this model, finite temperature effects 
lead to deconfinement.— In spin models of Mott insula- 
tors, a finite temperature deconfinement transition lib- 
erates spinon excitations. Recently, thermally induced 
spinon deconfinement was shown to occur in a class of 
pyrochlore antiferromagnets.— Entropic effects led to a 
finite temperature phase with algebraic correlations. 

We next derive the one-loop effective potential for the 
classical dimer model. Within the height representa- 
tion, this amounts to the effective potential for the two- 
dimensional sine-Gordon model. This can be readily 
acheived by shifting the field h and integrating out the 
Gaussian fluctuations. The effective potential obtained 
in this way is given by 



Ucsi^p) = -Ccos(27rA/K(^) 
cos (27rV^(p) In cos (27rV^(^) 



(38) 



where A is an ultraviolet cutoff proportional to the in- 
verse of the lattice spacing. In order to eliminate the 



dependence on the cutoff, we renormalize the above clas- 
sical effective potential by demanding that C7"g(0) equals 
Xj^ — Att'^kC,. Here, A^) is the correlation length as- 
sociated with the Debye-Hiickel approximation to (|37p . 
which consists in making the Gaussian approximation 
cos(27r-/K/i) « 1 — 27r^K/i^. Since 



1-1- 



■In 



(39) 



we have therefore that the renormalization condition im- 
plies \n{A-K'^ kC, / A?) = —1, such that the effective poten- 
tial becomes 



U^s{v) = ~C (1 - ^) cos (2^V^^) 

— cos {2n^/Kip) In [cos {2'IT^/KLp)] . 



This minimum energy density 



^0 = C7cff(0) - C(^k/2 - 1), 



(40) 



(41) 



changes sign at Kc = 2/7r, i.e., the usual critical value for 
the superfluid stiffness in the KT transition. 

That is, the VBS/plaquette system melts from an or- 
dered crystal to a critical phase with algebraic dimer cor- 
relations at the KT temperature 



KT 



(42) 



From Eq. (f38|) . we see that we can define an effective 
fugacity given by 



Ccff = C 



1 



(43) 



This immediately leads to RG equation for the dimen- 
sionless coupling C, — Ccff/A^, since up to terms of order 
higher than one loop, we can write 



A 



dA 



(44) 



so that 



(45) 



which is precisely the well known flow equation for the 
fugacity in the KT transition.*^ 

We now make contact with our relation of Eq. and 
a particular limit where its meaning becomes physically 
transparent. To this end, we consider the limit of large 
p that physically corresponds to a large negative v for 
which the columnar states are the ground states. In this 
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limit, the value of Tkt should tend to a divergent value 
in the limit p — > oo (or v — > — oo). Physically, for \v\ 3> t 
(i.e. , ignoring t) , the energy cost of breaking the columnar 
order in order to allow for a liquid state scales with 
itself. 

For |w| i and v < 0, there are 4 ground states — the 
columnar ground states. The corresponding phases break 
a discrete symmetry. Similar to the two-dimensional 
Ising modelr^ these discrete symmetries can and will 
be broken at finite temperatures. Domain wall energy 
penalties (scaling as \v\ (or p) multiplied by the perime- 
ter of the domain walls) will, at sufficiently low tempera- 
ture, overcome entropic contributions to the free energy 
and will force the system to order. Similarly, for any p 
(positive or negative), in the opposite regime of infinite 
temperature, the resultant classical random dimer cover- 
ing problem exhibits critical correlations.— The KT type 
transition derived above separates the critical high tem- 
perature and low temperature discrete broken symmetry 
regimes. If, unlike the QDM, the Hilbert space would 
be restricted to singlet dimers only in a low energy sec- 
tor and allow for unconstrained configurations at higher 
temperatures, then the critical phase found would ap- 
pear only in an intermediate temperature reg ime, such 
as may occur in spin systems, e.g., Ref. I 111 . In such 
an instance, for high temperatures, the correlations fall 
off exponentially. For temperatures in which the allowed 
configurations are restricted to dimer models, we will find 
algebraic decay at high temperatures (as in the infinite 
temperature limit of the QDM) and at yet lower tem- 
perature, the system would order into a VBS/plaquette 
phase which breaks discrete lattice symmetries. Such a 
behavior is found for Z^v model with large, but finite N. 
In such Zjsf models, a KT phase is found with both an 
upper and lower KT transitions The behavior of the 
zero temperature QDM, as we will elaborate on later, is 
that of — > oo. 



B. Finite temperature, p = 

The case of p = at finite temperature is less obvious. 
In this case, the finite temperature theory is given by 

/:classical = ^i^^^f - Cc0s(2^/l), (46) 

where A — K/T. The higher order derivative makes the 
upper and lower critical dimensions higher than in the 
standard sine-Gordon theory in Eq. ()37|) . which has an 
upper critical dimension =2. In the case of Eq. 
(|46)l the upper critical dimension is 13+ ~ 4, while the 
lower one is = 2. Therefore, provided the system is 
at the RK point and finite temperature, no phase tran- 
sition should occur as the finite temperature system lies 
effectively at the lower critical dimension, due to dimen- 
sional reduction. In order to see more concretely that 
no phase transition occurs in this case, let us consider 
the Renormalization Group flow equations in the general 



case of D dimensions instead of focusing solely on two di- 
mensions. In principle, a term ~ (V/i)^ is generated by 
fluctuations, but the inclusion of this term simply corre- 
sponds to a renormalization of p and we can still demand 
to be at the RK point by just requiring that the renor- 
malized p vanishes. The way of doing such a calculation 
has been discussed before4^ The same method can be 
used here. The resulting flow equations read 



dl 



C^ + iD- 4)A- 



dC 
dl 



^_,_(£^4)rCD^- 



c. 



(47) 



(48) 



We see that the second term between brackets in Eq. 
PS)) diverges \\i D — 2 dimensions. Therefore, the sys- 
tems does not reach any fixed point or line in this case. 
This more refined analysis clearly shows that no phase 
transition happens at a finite temperature RK point. We 
note in passing that for D = A the above flow equations 
lead to a four-dimensional KT transition. Such a tran- 
sition has been discussed some time ago in the context 
of quantum gravity;^ where a model similar to (|46p was 
considered. 



C. Zero temperature, arbitrary p 

For T = and p > there is also no phase transition. 
As we have already discussed, the higher gradient term 
becomes in this case RG irrelevant. It is well known 
that the model is in this case equivalent via duality to 
an electromagnetic (Maxwell) theory in 2 -I- 1 dimensions 
with a compact gauge field^S in which the charge-less, 
spin carrying, elementary excitations of the system, the 
so called spinous, are permanently confined. We have 
described above how spinon deconfinement can arise due 
to thermal effects in the gauge theory. 

Note that the absence of phase transition for p > Q 
and T = follows from an analysis of the effective height 
model. Thus, the effective model predicts a single VBS 
or plaquette phase (depending on the sign of z) for p > 
and vanishing temperature. However, it is known that 
for p > there is a first-order phase transition from the 
VBS state at larger p to a plaquette state at smaller p;^ 
see Fig. [TJ Up to now we have mentioned both phases for 
p > only in terms of the sign of z. In order to capture 
the phase transition, a numerical analysis of the lattice 
gauge theory of Eq. ([TB]) would be necessary. This will be 
done in a future publication. As a way of making progress 
by pure analytical means, we can alternatively consider 
higher harmonics of the cosine term in the Lagrangian 
([3T|) as derived from the action of Eq. ([30|) . This will lead 
to consider a Lagrangian that is discussed in Ref. [l^, 



C 



'-{drhf + 7: {^hf - zi cos(27r/i) - Z2 cos(47r/i) . (49) 
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If both zi and Z2 are positive, then as we discussed earher 
the higher cosine harmonic does not play any significant 
role at large distances, since it is RG irrelevant.— Thus, 
in such a situation we have just the plaquette phase as 
before when z > 0. However, if zi < then the term 
— Z2 cos(47rft,) is no longer irrelevant if Z2 > and this 
higher order term adduced from Eq. (f30|) must be kept. 

The application of mean-field theory to the Lagrangian 
of Eq.(l49)l leads to the solution 



ho 



1 

2^ 



(50) 



By computing the Gaussian fluctuations around /iq, we 
obtain a modified Debye-Hiickel approximation leading 
to a real-time height susceptibility of the form, 

1 

{h{uj,p)h{-uj,-p)) = - 



PP 



TT^Zn 



\4 



16z|) 



(51) 

Note that this type of Debye-Hiickel approximation is 
different from the one we have considered before, where 
the height field mean-field was an integer. For Z2 > 
and (zi/z2)^ > 16 we have a gapped spectrum for height 
field excitations. When Z2 > and (zi/z2)^ < 16, on the 
other hand, we have a pole for cj = at 



PZ2 



-(16z| 



zl). 



(52) 



These are two distinct regimes in a mixed phase between 
the plaquette and the VBS phases,^- which was recently 
shown to exist Another regime inside the mixed phase 
is the one where zxjzi — ±4, for which the excitations 
are gapless and = n/2, with n e Z, i.e., ho is either 
an integer or a half odd integer. The existence of gap- 
less exitations in the mixed state is important because 
it demonstrates that there is indeed a phase transition 
between the plaquette and VBS phase. It can now be 
seen that when Z2 — and zi = z a simple sign change 
in z does not lead to any phase transition between the 



plaquette and VBS phases. Indeed, after defining the 
couplings /C — {y/cp)~^ and y = ^/cpz, the RG equations 
would be in this case (the derivation can be easily done 
with the methods of Ref. 'i^ 



dy 
dl 



2/ 3 



dl 



= ic- 



(53) 



(54) 



Thus, irrespective of the sign of y (or z; note that dy/dl 
also changes sign), the RG equations remain invariant 
and no nontrivial fixed point can be found. The RG 
approach shows clearly that a sine-Gordon theory in 2 + 1 
dimensions does not exhibit any phase transition. This 
result is independent of the sign of z. This is not the case 
with the theory of Eq. digi) . 



Staggered 



Plaquette 



VBS 



p< RK point 



p> 



FIG. 1: (Color online) Schematic representation of the main 
different phases arising in two dimensions and at zero temper- 
ature. There is numerical evidence for a mixed phase between 
the VBS and the plaquette states (not shown here).— The re- 
sulting phase transition is discontinuous. 



Let us now turn back to the model with a single cosine 
harmonic and ^ 0, Eq. (|3ip . at zero temperature. 

By integrating out the Gaussian fluctuations in the height 
field 



K{c + pp^) +p'^ + Mfj cos{2-K^Kip) 
K{c ijj^ -I- pp^) +p'^ 



(55) 



where similarly to the finite temperature case, — An'^Kz gives the inverse square of the Debye-Hiickel length. We 
obtain 



'Z ( 1 - -\ — I cos(27r\/X^) - — J - cos{2TT\/~Kip) 
\ 8 y c I 16 yc 



Ik Kp^ - 167r2zcos(27rViir^) 
64^? 



In 



Kp + An^/ Kz cos(27r\/^<^) 

1A2 



(56) 
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By using the renormalization condition [/"fj(0) = we finally obtain the renormalized effective potential 

fW Kp"^ - IQt:'^zcos{2t:^/Klp) 



-zil - -\ — I cos{2T:\fKLp) - — J - cos{2TT\/~Kip) 
\ 8 y c I 16 yc 



647r 



■In 



Kp + Any Kz cos(27rv K^p) 
Kp + At^^TKz 



(57) 



At the RK point we have p = and the above effective 
potential becomes 



f/eff(<y3) 



-z \ 1 




cos(27rv Ki^') 



— \\ — cos(27r-\/i?95) In cosil-K^fKLp) 
° " c 



(58) 



Interestingly, Eq. ([55)1 gives the same type of effective 
potential as in Eq. (|40)) . which was obtained for p > 
and high temperatures. Note, however, the presence of 
the extra parameter c in Eq. (|55|) . On the other hand, 
if we introduce the dimensionless parameters c = c/A 
and K = K/A along with (j) = V^if, we see that for the 
special choice c—l/K the effective potential becomes 



■kK 



Ues{(f>) = 1 " ^ cos(27r 



TTZK 



cos(27rVif0)ln cos(27rVif0) 



(59) 



which has now exactly the same form as (l40l) . The above 
effective potential corresponds to a KT-like transition 
with a critical stiffness Kc — S/tt. Indeed, the minimum 
energy density is given by 



Eo = i7eff(0) 



-z 1 



irK 



(60) 



which changes sign precisely at Kc — 8/tt. A KT-like 
transition at zero temperature was also obtained recently 
for a QDM.^i^ 

There is a deeper argument behind the choice c = 1 /K. 
It actually corresponds to a self-dual point of the model 
in a special limit. To see this, first we note that for p = 
the Lagrangian ((3T|) is dual to the U(l) symmetric lattice 
Lagrangian^ 



C/(l) 



E 



2 2c 



2nmj 



(61) 



where ni,mi € 1j and —it < 9i <tt. The self-duality and 
its relation to a KT transition in 2 -f- 1 dimensions follows 
from the dualization of the following Zjv model. 



C2 



E 



K /27r 



1 



2 

27r 



N 



2a U^ '^^ 



27rm, 



(62) 



in the iV ^ oo limit, where qi — 0, 1, 2, . . . , TV — 1. In 
the limit iV — > oo the Lagrangian ([62]) becomes ([6T|) . As 
shown by Amit et al.p^ a duality transformation relates 
the partition function of the model to its dual as 



Z(k, l/c) - Z{N^c/{A7T^), N^/{4Tr^K)), 



(63) 



and we see that the couplings are not only inverted by 
duality, but also exchanged. Therefore, we obtain the 
self-duality relation c/K = Att'^/N^. For large N the 
self-dual point with c — 1 /K leads to a phase transition 
as K is varied. This self-duality point is precisely what 
we have demanded to obtain Eq. ([5^ . with a resulting 
KT transition. The self-duality breaks down along the 
line c = K . In this case no phase transition occurs. 

Note that in spite of the similarities between the de- 
scribed self-dual quantum critical point and the usual 
KT transition at finite temperature, there are actually 
important differences, making them distinct universality 
classes. That is the reason why we refer to the self-dual 
quantum critical point as "KT-like" rather than simply 
"KT". In order to emphasize this point, let us consider 
the spin-wave theory of the the model (|6ip in the con- 
tinuum limit and at the self-dual point. We have the 
correlation function between the XY spins: 

C(r,r) = l^^^e{r.r)^-^eim) 

= AAexp|-i[G(r,r)-G(0,0)]|, (64) 
where A/" is a normalization constant and 



(65) 



In the static limit we have a regime looking very much 
similar to KT, since 
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G(0,r) 



An V4a/ 



const, 



(66) 



where a is a cutoff. However, tliis leads to a different 
anomalous dimension for the correlation function C(r, r) 
in the limit r ^ 0, i.e., 



V 



AnKr 



1 

32' 



(67) 



which is to be compared to ry = 1/4 for the usual KT 
transition. 

On the other hand, note that for large r the Green 
function behaves like 



G(T,r) 



■ exp 



Attt '"^ \ 2t2 
which can be derived from the exact representation 



(68) 



G(r,r) 



1 r du 



ur 



2(r2 + Av?) 



X exp 



2(t2 + 4u2) 



(69) 



where /q is a modified Bessel function of the first kind. 
For large distances with finite r we have an asymptoti- 
cally static behavior: 



(a) z > 



Plaquetle 



KT critical point 



(b) z < 




K 



FIG. 2: (Color online) Schematic representation of the phases 
of the two-dimensional QDM at the RK point. The zero tem- 
perature KT (quantum) critical point is given by Kc = S/tt. 
Note that for z > [Panel (a)], the KT-like transition is 
from a plaquette to an RVB state, which is a liquid of valence 
bonds. In Panel (b) the RVB state emerges out of a VBS 
state at 2 < 0. 



KT transition line 




Staggered 



Plaquette 



RK poitn 



G{t,v) 



(70) 



The zero temperature transition of the two- 
dimensional QDM describes a phase transition where 
a plaquette or VBS state changes into a spin liquid, 
here thought as a liquid of valence bonds. A very 
schematic representation of the phase transition at 
zero temperature is shown in Fig. [51 As in the finite 
temperature KT transition of the QDM, we have here 
that a VBS state corresponds to a "superfluid" phase, 
while the spin liquid will be the "normal fluid" of the 
transition. 

Putting all of the pieces together, we arrive at Table 
I that encapsulated our results at the beginning of this 
section. We obtained these results from the analysis of 
effective height models via duality arguments and the 
study of effective potentials. A schematic phase diagram 
in the pT-plane is shown in Fig. [3) 



FIG. 3; (Color online) Schematic phase diagram of the two- 
dimensional QDM in the pT-plane. For large p, the KT 
temperature is high and varies linearly with p with a slope 
of (2/7r) (Ea. (l42|) '). Here we have also indicated the mixed 
phase.— Temperature effects melt the dimerized phases to 
a valence-bond (VB) liquid phase, which is actually the same 
as a thermally induced RVB state. 



arises via duality. The effective model in this case is a 



compact gauge theory 



2.32 



C = ^-{drA-f + ^(V X Af + i^(V X V X A)^ (71) 

As before, the RK point corresponds to p = 0. We now 
derive the corresponding dual model at the RK point. 
This is easily done by considering the lattice model ver- 
sion of (j71[) in the Villain approximation. In this approx- 
imation, the lattice action reads 



VI. THREE-DIMENSIONAL EFFECTIVE 
MODEL FOR THE RK POINT AND DUALITY 



S = 



Now we will consider the (3 -f l)-dimensional quantum 
dimer model. In this case no pure scalar field theory 



(72) 
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where nj and xiij are integer fields and the lattice sites 
j are in 3 + 1 dimensions. The first important step is 
to introduce new integer-valued auxiliary fields using the 
Poisson formula. This leads to the action 



S 



E 



V X V X Aj] 



(73) 



where Mj and Nj are integer fields. As the gauge field 
Aj is compact, integrating it out yields the local con- 
straint 



V^Mj + V^Nj- - V(V • Nj) = 0. (74) 
This constraint is solved by the equations 



Mj = V X V X Lj, 



(75) 



Invoking the Poisson formula once more to promote the 
integer field Lj to a real field a.j , we obtain the action of 
the dual model in the form 



S 



E 



K 



(Vraj 



1 



(V X V X aj)^ — 27rmj • a.j 



(76) 

where Wj is an integer-valued field. Two important re- 
marks are in order. First, we note that, as before, the 
couplings are inverted and exchanged. However, the 
gauge field appearing in the dual model is not compact. 
Second, there is a coupling between the gauge field and 
an integer field. To understand the physical meaning of 
this coupling, note that due to gauge invariance in space 
direction^'* the constraint V • rij = holds. This means 
that we can sum over field configurations with loops (akin 
to vortex loops) j^i^iS This constraint enters in the par- 
tition function as a Kronecker delta. Employing the in- 
tegral representation of the Kronecker delta and partial 
summation yields 



E 



Y(V,a,)2 + l(Vx Vxa,f 



+ ^(V0J - 27raj-) • nj] , (77) 

where Oj g [0, 27r]. In the sum over rij, we introduce in 
the action a soft vortex-core term ecn^/2 with a small 
core energy ec-^^ By the Poisson summation formula, 

^ = E 



Y(V,a,)2 + l(Vx Vxa,)2 



where ■0 is a complex scalar field. Quantum fiuctuations 
generate a term ^ ln(^A)(V x a)^ (with a positive coeffi- 
cient), where ^ is the correlation length. This term ren- 
ders the double-curl term irrelevant at large distances. 
Quantum fluctuations also generate a |9r'0P term. The 
resulting effective action is similar to that of an Abclian 
Higgs model in four spacetime dimensions. The rcnor- 
malization group (RG) flow for such a theory is known^ 
to exhibit a Gaussian fixed point with a runaway flow. 
The same occurs here and no nontrivial fixed point is 
found. Indeed, the coupling is marginally irrelevant. 
This leads to a first-order phase transition. Thus, no 
quantum critical point appears in this case. 

An interesting aspect of this first-order phase transi- 
tion is that the mass gap of the gauge field is induced 
by fluctuations, through the so-called Coleman- Weinberg 
mechanism i^Siffi Here we have an interesting modified 
version of it. This can be seen more straightforwardly 
by integrating out the gauge field while neglecting the 
scalar field fiuctuations, i.e., assuming ip uniform. This 
yields the effective potential. 



(80) 



where A is an ultraviolet cutoff. Note that, in contrast 
with Ref. [5^, there is a prefactor of j^/'P instead of a pref- 
actor I '01'^ that multiplies the logarithmic term. This is a 
consequence of the p* term in the gauge field propagator. 

The first-order phase transition disorders a three- 
dimensional VBS into a RVB state. Recall that in two 
spatial dimensions the quantum phase transition from a 
VBS to a RVB state exhibited a quantum critical point 
governed by a KT transition. 

It is also useful to consider the case p > and K ^ 00, 
in which case we have instead Eq. ([72ll the lattice action. 



^ = E [^i^rA, - 2nn,f + ^{\/ x A, - 2^m,f 

The first step of the duality transformation leads in this 
case to an action featuring integer fields M^- and Nj of 
the form 



-JE^m, 



J p J 



(82) 



-t- — {VOj - 2TTSLj - 27rmj)^ 



(78) 



In the continuum limit, the dual model (|78p reads 
C = :|(a,a)2 + l(Vx Vxa)2 



where the constraint 



VrM..j = V X Nj, 



(83) 



+ \{V-2^^^)^P\'+r\i^\' + -\^\\ (79) we obtain finally 



must be fulfilled. The constraint is easily solved by in- 
troducing a new integer field such that Mj = V x Lj 
and N, = VrL,-. With the help of the Poisson formula 



13 



£(Vxa,)^ 



(84) 

Here the field is a real field while rij is an integer field 
and we have performed a trivial rescaling. The above 
action is similar to the action (|76p . except that we have 
here a normal curl term instead of a double-curl one. 
Therefore, the continuum limit of the dual lattice model 
above reads 



C 



1 



p>0 



P_ 

2c 



(V X a)2 



\{V-2tii^P a)V|' 



r\i^\^ + ^\^t. (85) 



At zero temperature the above Lagrangian describes ef- 
fectively the same first-order phase transition as in the 
case of the theory (l79l) featuring a double-curl. There- 
fore, we see that for p > no quantum critical point 
exists. The first-order phase transition creates a three- 
dimensional liquid of valence bonds out of a VBS state. 

At high temperatures, on the other hand, the situation 
is very different, since in this case a(T, x) « a(x). and 
the term [drSiY can be neglected. By defining new fields 
according to = ^jjy/T and h = ayj pl{cT), we obtain 

1 V 

^P>o « -(Vxh)2 + |(V-zeh)vI/p+r|vI/|2 + -|vI;|4, (86) 



where e = 2n\' cT and v = uT . The above expression 
is exactly the same as the Ginzburg-Landau (GL) free 
energy of a superconductor in three dimensions. In this 
case, it is known that the order of the phase transition 
is determined by the ratio w = v/{2e^)§ ^-^^ For w > 
0.55, the phase transition is second-order jSiSi otherwise 
a first-order transition takes place. Duality arguments 
and Monte Carlo simulations^^'^^ indicate that when w > 
0.55, the correlation length and specific heat exponents 
attain XY values, otherwise the phase transition is a first- 
order one. 

At high temperatures a similar analysis can be made 
in the dual model ((79)) to obtain the Lagrangian 

i(VxVxh)2 + |(V-ieh)1'|2+r|*|2-f-^|5'|4, (87) 



where, as before, e = 2'KycT and v = uT . If it was 
not for the double-curl term, the above expression would 
be exactly the same as the free energy for a supercon- 
ductor in three dimensions. However, we can proceed 
similarly as in the zero temperature case and argue that 
thermal fluctuations will generate a (V x h)^ term with 
a positive coefficient, making once more the double-curl 
term irrelevant. Note that this time the generated curl 
term does not contain any logarithm in front of it, since 
the high temperature theory is three-dimensional. Thus, 
by taking thermal fiuctuations into account, the finite 
temperature theory behaves in the same way as a GL 



superconductor in three dimensions. This result approx- 
imately agrees with recent simulations of a classical dimer 
model.— However, recent results suggest a different uni- 
versality class.--"-'* In fact, our argument of irrelevance 
of the double-curl has to be taken with care. In order to 
appreciate the subtleties of the problem, let us consider 
a simple one-loop calculation. In this case, the gener- 
ated single-curl term is ~ Ce^(V x h)^, with a correlation 
length ^ ^ (r— Tc)"'', where Tc is the value of r at the crit- 
ical point and v is the corresponding critical exponent. 
Thus, it seems that for fixed and very near to r^, the 
coefficient of the generated (V x h)^ term is large. There- 
fore, for large distances (corresponding to |p| — > 0) the 
double-curl term seems to be indeed irrelevant. However, 
this argument may be too simplistic for the following rea- 
son. Simple dimensional analysis of the Lagrangian (j87p 
shows that has dimension of (length) Thus, its 
renormalized counterpart should behave near the critical 
point as e|j ~ so that the actual one-loop renormal- 
ized coefficient of the fluctuation-generated (Vxh)^ term 
behaves, after replacing by e^, as ~ Therefore, 
precisely at the critical point this term vanishes and the 
double-curl term plays again a role. However, if we start 
with the theory already at the critical point, i.e., with 
r = 7'c, it is easy to determine the anomalous dimension 
r\h of h to all orders in peturbation theory, provided a 
perturbative fixed point is found. The method resem- 
bles the one used in the study of critical fluctuations in 
superconductors.^^ In the present case it yields r\h = 3, 
such that the gauge field propagator has the behavior 



(^i(p)^j(-p)) 



1 



1 



which is identical to the infrared behavior of the GL 
modeli^ This highlights, once more, the irrelevance of 
the the double-curl term at large distances. Thus, we see 
that there are two possible situations relative the critical 
point, that correspond to two different fixed points. This 
rough analysis shows that the double-curl term in Eq. 
(|87|) is dangerously irrelevant. Thus, in order to unam- 
biguously determine the universality class of the model of 
Ea. (j87p . a more quantitative analysis is necessary. This 
will be subject of a future publication. 

While in the case of the Lagrangian of Eq. ([87]) the 
dangerously irrelevance of the double-term makes it ques- 
tionable to discard it at the critical point, the same prob- 
lem does not arise with Eq. ((7T|) for a fixed p > 0. In 
that case the double-curl term can safely be neglected at 
large distances and the lattice model (I8ip corresponding 
to X ^ oo can be used. 
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VII. CONCLUSION 

In this paper we have derived an extended lattice gauge 
theory action for the QDM. From it we were able to de- 
rive effective models to study the phase structure of two- 
and three-dimensional QDMs. For the two-dimensional 
case we have calculated the one-loop effective potential. 
At the RK point the effective potential is characteristic of 
one exhibiting a line of fixed points leading to a KT-like 
transition. This KT-like transition in 2 -I- 1 dimensions 
and zero temperature is a consequence of the anisotropy 
and self-duality of the effective model. It governs a uni- 
versality class which is actually distinct from the tradi- 
tional KT transition. This is also in stark contrast with 
the (2 -I- l)-dimensional KT-like transition discussed in 
the past in the context of relativistic U{1) spin liquidsj^ 
In that case the system was isotropic and the propagator 
non-analytic. It is still controversial whether a (2 -I- 1)- 
dimensional KT transition may occur in this casci^ 

The quantum critical point of the two-dimensional 
QDM was related to a self-dual point of a family of 
theories in the limit N oo. This self-duality is strongly 
dependent of being at the RK point. One interesting as- 
pect in the duality scenario of the two-dimensional QDM 
is that the RK point is protected by the duality sym- 
metry. In other words, as a consequence of the intrinsic 
anisotropy, the parameter p ex. {t — v) does not get in- 
verted by the duality transformation, so that the RK 
point (i.e., the point for which p = 0) is preserved by the 
duality transformation. 

It is well known^ that in 1 + 1 dimensions the stan- 
dard Zn model (i.e., without higher gradients) undergoes 
a KT transition for N large. This is not surprising, since 
the Zn model becomes ioT N ^ oo the XY model. In 
our case, the KT transition occurs one dimension higher 
as a consequence of the higher gradient term and the 
anisotropy. To this KT transition there is also a corre- 
sponding Zn model, Eq. ([5^ . whose large A'^ critical 



point governs the KT transition of the height model at 
the RK point. As we have discussed, in this case there is 
an interesting additional feature, which is the self-duality 
at large A^. 

At finite temperature and above the RK point (t > v) 
a genuine KT transition occurs. In contrast with the 
zero temperature KT-like transition, the finite temper- 
ature transition is a consequence of dimensional reduc- 
tion induced by temperature rather than of the intrinsic 
anisotropy of the system. At the RK point, on the other 
hand, no phase transition occurs at finite temperature. 
In this case dangerous infrared singularities prevent a 
phase transition to happen. 

For the three-dimensional case we have derived the 
dual model at the RK point and shown that at zero tem- 
perature a first-order phase transition takes place. A 
second-order phase transition can only occur in this case 
at finite temperature. The finite temperature phase di- 
agram features first- and second-order phase transitions 
separated by a tricritical point, in a scenario very remi- 
niscent from phase transitions in superconductorsj ^^'^^i^^ 

Finally, it should be noted that the approach employed 
in this paper can be used to derive far more accurate re- 
sults than those that we have discussed. In particular, 
the phase boundaries near the RK point in the schematic 
phase diagram can in principle be determined more pre- 
cisely via Monte Carlo simulations of the extended lat- 
tice gauge theory introduced in Sect. III. This will be 
the subject of a future publication. 
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